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ABSTRACT 


A  computer  program  to  calculate  the  mass,  center  of  gravity  location, 
and  moments  of  inertia  of  a  system  of  coaxial  bodies  of  revolution  is 
presented.  The  derivation  of  equations  used  by  the  program,  instructions 
for  setting  up  inputs,  and  a  sample  case  are  also  given. 
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LIST  OF  SYMBOLS 


a0*  al 


y-xntercept  and  slope  of  a  straight  line,  respectively 
transverse  moment  of  inertia  about  the  center  of  gravity 


rx>  Iy.  I2  moments  of  inertia  about  the  x,  y,  and  z  a-:es 


polar  coordinates  used  in  transformation  of  y-z  plane 


x,  y,  2 


V  xf 


Subscripts 


radius  of  a  circular  arc 

right-handed,  orthogonal  coordinate  system 

center  of  gravity  position  along  x-axis 

lower  and  upper  bounds,  respectively,  of  surface  along  x-axis 
density 


coordinates  of  center  of  circular  arc 

value  for  segment  of  body 

total  value  for  body 

upper  and  lower  surfaces,  respectively 
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I .  INTRODUCTION 


When  designing  a  projectile,  one  must  consider  not  only  the  exterior 
configuration  of  the  body,  but  its  physical  properties  as  well,  since 
these  will  directly  influence  the  flight  behavior  of  the  shell.  By 
physical  properties  we  mean  mass,  center  of  gravity  location  and  the 
axial  and  transverse  moments  of  inertia.  It  is  possible  to  compute  these 
properties  manually,  but  this  task  for  a  relatively  complex  projectile 
is  a  very  tedious  one.  The  program  described  in  this  report  (coded  by 
D.  Solmon)  enables  the  designer  to  obtain  accurate  values  for  the  physi¬ 
cal  characteristics  of  his  designs  with  minimum  effort. 

Minimization  of  user  effort  necessarily  implies  some  constraints. 
However,  the  constraints  to  be  applied  must  not  seriously  degrade  the 
ability  of  the  program  to  handle  complex  bodies.  With  this  in  mind,  one 
major  assumption  was  made  in  designing  this  program;  namely,  that  objects 
to  be  considered  by  this  program  will  consist  of  coaxial  bodies  of 
revolution  only.  More  generalized  programs  are  available  which  handle 

the  asymmetric  case,  but  which  also  require  more  work  on  the  part  of 

1* 

the  user 

This  report  presents  the  derivation  of  the  equations  used  by  the 
program  and  instructions  for  setting  up  the  inputs.  A  sample  case  is 
included  for  illustrative  purposes.  A  complete  listing  of  the  program 
with  all  subroutines  may  be  found  in  the  Appendix. 

II.  DERIVATION  OF  EQUATIONS 

Consider  an  axisymmetric  shell  of  uniform  density  (Figure  1)  having 
the  x-axis  as  its  axis  of  symmetry.  The  shell  is  bounded  radially  by  rf 
and  r  (where  r  and  r  are  functions  of  x  and  r„  =  y.  and  r  h  y  in 
the  x-y  plane)  and  in  the  x-direction  by  xq  and  x^  (where  xo  is  not 
necessarily  located  at  the  origin  as  shown  in  Figure  1). 


* References  are  listed  on  page  21. 
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A.  Mass  and  Center  of  Gravity  Location 


In  general,  we  know  that  for  uniform  density 
m  =  /  dm  =  p  /  /  /  dx  dy  dz. 

Transforming  the  y-z  plane  to  polar  coordinates  and  setting  the 
limits  of  integration,  we  get 

dy  dz  =  r  dr  d6. 


which  reduces  to 


Ke  also  know  that 


m  =  p  f 


f  / 2*  /  u  r 


x  o 
0 


dr  dB  dx, 


m  =  ’p  /  (y  -  yr  )  dx, 
x  u  1 
o 


/  x  dm 


eg  /dm  ' 


therefore. 


/  (>•  2  -  y„2) 


/  f  (y  2  -  y  2)  dx 
x  v/u 
o 


To  calculate  the  total  mass  and  center  of  gravity  location  for  a 
composite  body,  we  first  calculate  the  mass  and  center  of  gravity  loca 
tion  for  eac^  section  (m,  x  )  and  use  the  following  relations: 


R.  Moments  of  Inertia 


The  principal  moments  of  inertia  can  be  defined  by 

I  =  p  /  /  f  (y2  +  22)  dx  dy  d2, 

I  =  p  /  /  /  (x'  +  2 2 )  dx  dy  dz, 

I..  -  p  /  S  f  (x2  +  y2)  dx  dy  dz, 

z 

where  1,1,  and  I  arc  the  moments  of  inertia  about  the  x,  y,  and  z 
x  >’  2 

axes  respectively,  Tor  an  ax i symmetric  body  whose  axis  of  symmetry  is 
the  x-axis, 

1=1. 

y  = 

1 .  Axial  Moment  of  Inertia  . 

I  =  p  /  /  j  (>'■’  +  z2)  dx  dy  dz. 

Transforming  to  polar  coordinates  and  setting  limits  of  integration, 

x ..  ....  >• 

I  =  /  /'  /  U  r3  dr  dO  dx 

X  X  C  V  . 

O  '  .. 


there foi  c , 


nr 


Xf 


I  .  -  -  >  r  4 )  dx, 


for  a  com])ositc  body. 


"  I 

i 


r  x . 


Transverse  Moment  of  Inertia. 


therefore, 

2  I  =  I  .  +  I_  =  ■  j  f  (2x-;  +  y?  +  z')  dx  dy  dz. 

I  =  p  f  f  f  x 2  dx  dy  Jz  +  f  f  f  (y2  +  z2)  dx  dy  dz 

=  ,  if  /  x2  dx  dy  dz  +  4  I  . 


Transforming  to  polar  coordinates  and  setting  limits  of  integration. 
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ly  ~  J  !x  +  P  fj  ir  X1  T  dr  de  dX» 


x  o  y. 
o  ’  a 


therefore, 


1 


I  =  t  I  +  no  /  ^  (y  2  -  y  2)  x2  dx. 

\/  /  Y  ^  '••'ll  '  0  ' 

o 


y  2  x  x_  v'u  'H 


For  a  composite  body,  the  total  moment  of  inertia  about  the  y  (or  2) 
axis  is  given  by 

>y  *  I  v- 

't  i  yi 

To  get  the  total  transverse  moment  of  inertia  about  the  center  of 
gravity,  we  use 


I  =  I  -  nr  x 


ch  yt 


t  c8» 


III.  USH  OF  THC  PROGRAM 

As  previously  stated,  this  program  assumes  a  system  of  coaxial  bodies 
or  "evolution.  If  the  x-axis  is  chosen  as  the  axis  of  symmetry,  then  the 
surfaces  of  the  body  may  be  generated  by  rotating  y  =  f(x)  about  the 
x  xis  The  origin  is  usually  taken  (a)  at  the  nose,  with  the  positive 
x-a^cs  pointing  rearward  or  (b)  at  the  base,  with  the  positive  x-axis 
pointing  forward.  The  center  of  gravity  is  computed  from  the  chosen 
origin.  As  presently  constructed,  the  program  will  handle  two  types  of 
functions;  circular  arcs  and  straight  lines.  Circular  arcs  are  of  the 
form 


y  =  yc  f  '  r2  -  (x  -  xc)2, 

where  (xc>  y  )  is  the  location  of  the  center,  and  R  is  the  radius.  Some 
care  must  be  taken  to  insure  that  the  quantity  under  the  radical  sign  is 
never  negative  in  the  applicable  x-interval .  (The  quantity  could  go 
negative  near  |x-xcJ  =  R,  due  to  round-off  errors.)  Taking  the  origin  at 
the  nose  will  usually  circumvent  this  problem. 
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Straight  lines  are  of  the  form 


y  =  a  +  a  x, 

•A.  1 

'w  i 

where  and  a^  are  the  y-intcrccpt  and  slope  respectively.  Associated 

with  each  function  i  •>  an  interval,  (x^,  x^.1  ,  within  which  it  is  applicable 

and  the  density,  p,  of  the  area  lying  immediately  below  the  function  with¬ 
in  the  interval. 

bach  function  is  input  on  a  single  data  card.  The  cards  may  be 
arranged  in  any  order,  with  a  blank  card  fol lowing  the  last  data  card 
ot  a  case.  Cases  may  be  stacked.  The  data  cards  are  of  the  following 
tori::,  with  dat.i  1  i  e  1  d  --  '■>«.'  i  ;i  ;■  f  <  ; .  coh.sin-.  Iv  r:.'i  1  point  mo  T  ' 1 : 1 1 .  •  i ,  ■  •-. !  . 

A.  ci  reel  a  r  rc  s 


blank  for  y 


K.  Straight  Lines 


Card 

Columns 

1-10 

11-20 

21-30 

31-40 

41-50 

51-60 

61-78 

79 

80 


Content 


blank 


P 

alphanumeric  code  for  identification  of  output 

blank 

2 


Care  should  be  taken  to  insure  that  the  units  of  measurement  used 
for  density  are  consistent  with  the  units  of  length  used  on  the  drawing 
from  which  the  functions  were  derived. 


The  program  prints  out  the  input  data,  for  checking  purposes,  as  well 
as  the  computed  values  of  mass,  center  of  gravity  location,  and  axial  and 
transverse  moments  of  inertia.  The  units  of  these  computed  values  are 
dependent  upon  the  units  of  the  input  data. 


IV.  SAMPLE  CASE 

The  sample  case,  shown  in  Figure  2,  is  the  105mm,  HE,  Ml  artillery 
projectile  with  M73  dummy  fuze.  The  shape  is  rotationally  symmetric 
except  for  two  fuze  wrench  slots  on  the  M73  which  were  ignored  for  present 
purposes.  The  densities  of  the  various  materials  which  make  up  the  round 
are  listed  in  Table  1. 
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Table  I.  Densities  of  Materials  Used  for  Sample  Case 


Section 

Material 

Density 

(lb/in3) 

Fuze 

Steel 

.2833 

Body 

Steel 

.2833 

Base  Plate 

Steel 

.2833 

Rotating  Band 

Gilding  Metal 

.3128 

H.E.  Filler 

Comp  B 

.0549 

Tabulations  of  input  and  output  are  shown  in  Tables  II  and  III 
respectively.  A  comparison  of  computed  values  with  standard  values  is 
given  in  Table  IV.  The  computed  values  are  in  good  agreement  with  the 
standard  values  with  the  maximum  error  (*2.8%)  occurring  in  the  trans¬ 
verse  moment  of  inertia  computation.  Keep  in  mind,  however,  that  these 
standard  values  are  the  mean  values  of  measurements  taken  on  a  sample 
of  production  rounds  whose  actual  shapes  may  vary  slightly  from 
standard.  This  sample  case  was  selected  to  give  the  reader  an  idea  of  the 
degree  of  complexity  which  can  be  handled  by  the  program.  For  known 
shapes  with  known  densities,  the  computation  is  nearly  exact  (within  the 
tolerance  imposed  on  the  integration  routine). 
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Table  IV.  Comparison 

of  Standard 

and  Computed  Values  for 

Sample  Case 

Quantity 

Units 

Stand . 

Computed 

Error 

Co) 

Mass 

lb 

33.0 

33.087 

+0.3 

Center  of  Gravity 
(from  nose) 

in 

12.264 

12.291 

+0.2 

Axia]  Mom. 

lb-in^ 

79.488 

80.453 

+  1.2 

Trans.  Mom. 

lb- i n? 

770.803 

792.18 

+  2.8 
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C 

C 

C 


DIMENSION  IC  00E ( 40 )  i  XCI40),  YCUO),  RCI40).  A0140).  A1I4C), 

1  X0(ft0!,  XF(80I.  RCI  40)  ,  X0P(  80 1  t  X<80)»  XT<8C),  FYI401,  MNI40I, 

2  NFUNI40I.  XLI40I,  XU<40) 

COMMON  XC,  YC,  RC.  AO,  AL,  NuP,  NLO ,  FLIP,  ICCCE 
EXTERNAL  FX 
P I *3. 1 A  1 592654 

INITIALIZATIONS 

1  2F*0. 

XM«0. 

AM  I *0. 

ee*o. 

A0( l I =o. 

Aid  I  -0, 

XCI1I =0. 

veil) -0. 

RCI 1 ) =0. 

R0(  l I =0, 

XOI 1 t =0. 

XF ( 1  I =0. 

RE AO  INPUT,  REARRANGE  |F  NECESSARY,  AND  PRINT  OUT. 

00  5  1=2,100 

READ  (5,28)  XC ( I) ,YC(  I  l.KC  M  ), X0( I  I  ,XF  II  ) .ROII) , ! ACOOE(K) , K*  1,2), 
l  ICOCEMI 

if  i  icoocin.ec. o)  go  to  6 

A ACOC E * AC COE  I  1  I 

A8C0CE  =  ACCDE ( 2  I 

IE  (  ICOOEIU.NE.2)  GO  TO  2 

AO  I  I  I  =  XC  (  II 

Aim  *YC  I  I  > 

XF (I ) =X0 I  I  ) 

XOI I  I  =RC  I  I  ) 

xen  i  =0. 

YCII 1=0, 

RCI I) =0. 

CO  TO  3 

2  ao<  i  no. 
aii  i  no. 

3  IF  I  I , G T . 2 )  GO  TO  A 
WRITE  16,311 

A  WRITE  (6,321  A0(  I  I , A1  (  I  I ,XC<  I  I , YCC  I  1 ,RC!II , XOI  I  1 , XF  I  I  )  ,R0  (II  , 

1  AACCOC .ABCOOE 

5  CONTINUE 
N  =  20 

GO  TO  7 

6  M  =  I  -  1 

7  WRITE  (6,29) 

L  I=N 

ACCOE (  1  ) =AACOCE 
ACODE  <  2  )  =  A  OC  OCE 

DIVIDE  BODY  INTO  REGIONS  USING  BREAK  POINTS. 

00  8  I  =  1  ,  N 
X  (  I  )  =  X  0  (  I  I 
J  -M  I 

8  X ( J)  =  XF(  1  1 


HA  IN  1 
MAIN  2 
MAIN  3 
HA  IN  4 
HAIN  5 
HA  IN  6 
HAIN  7 
MAIN  8 
HA  IN  9 
HAIN  10 
HAIN  11 
HAIN  12 
HAIN  13 
HAIN  14 
HAIN  IS 
HAIN  16 
HAIN  17 
HA  IN  18 
HAIN  19 
HAIN  20 
MAIN  21 
MAIN  22 
MAIN  23 
MAIN  24 
MAIN  25 
MAIN  2 6R 
HAIN  27R 
HAIN  28 
MAIN  29 
MAIN  30 
MAIN  31 
HAIN  32 
HAIN  33 
MAIN  34 
MAIN  35 
MAIN  36 
MAIN  37 
MAIN  38 
MAIN  39 
MAIN  40 
MAIN  41 
MAIN  42 
MAIN  4  3  W 
MAIN  44W 
HAIN  45W 
MAIN  46 
MAIN  47 
MAIN  48 
MAIN  49 
MAIN  SOW 
MAIN  51 
HAIN  52 
MAIN  53 
HAIN  54 
MAIN  55 
MAIN  56 
MAIN  57 
MAIN  58 
MAIN  59 
HAIN  60 


( 

y 


r 
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onn 


NN*2  *N 

N«NN 

L*1 

00  12  1*1. N 
PLOP  *  0. 

CALL  PMIN  (X.NN.XMIN. ITI 
X8P(  L I »  X  M I  N 
IF  (  l.EO.l 1  GO  TO  9 
IF  (XeP(L).EO.XBPCL-m  PL  OP  *  1«  0 
9  J  *0 

00  10  K.l.NN 

IF  1 IT.EC.K)  GO  TO  10 

J*  J*  l 

X  (  J I  *  X  (  K  I 

10  CONTINUE 
NN«NN-l 

IF  (FLOP. 60. t.OI  GO  TO  11 
L*l+1 

11  CONTINUE 

IF  (NN.EQ.1I  GO  TO  13 

12  CONTINUE 

13  CONT  INUE 

IF  (X(l).EO. XBPlL-t  I  »  GO  TO  14 
XBPI L  )  *X (  1 ) 

GO  TO  15 

14  L  =L“  1 

15  CONT INUE 
X0(1 ) 1 X  BP ( 1  I 
X  F ( 1  ) =X«P(L I 
IM 

16  11*1 

x T  (  I  I * ( XBP (  I  I ♦ XB P (  I  *  l  I  1/2. 

X  *  1 

SEPARATE  REGIONS  INTO  LAVERS  OF  UNIFORM  DENSITY. 

CO  19  J* I .LI 

IF  ( XOi JI.GT.XTt  I  I.OR.XTI I  I.GT.XF I JII  GO  TO  19 
IF  (RC( JI.EO.O.  )  CO  TO  17 

F  Y(K  l*YC  (JI»SO«T  («C  I  J  l**2-l  XI(  I  )-XC(  J1)*»2MAC(  JI  +  A1  (  JI»X  T(  I  I 
IF  (  ICOOEI JI.NE.-I  )  GO  TO  18 

FY(X»*YC(Jl-S0Rr(PC(JI**2-tXT(I»-XC(J))*»2)»A0(J)*Al(J)*XKII 
GO  TO  18 

17  F  Y ( K ) *A0 ( J»*A1 1 J  )*XT  t  I  > 

IB  CONTINUE 
MNIKlsJ 
K*K.  1 

19  CONTINUE 
X  =  K-1 

J  J  =  1 

20  CONTINUE 

IF  IK, £0.21  GO  TO  22 
CALL  PMAX  (FY.K.XMAX. ITT 
Nf UN (  J J I =  “N(  I  T  | 

ICODEI  JJ)*IC0C£<  in 
J  =  0 

DC  ?1  I  x  =  1 , K 

IF  (  I  T. EC. IX  I  GO  TO  21 

J*  J»  1 

I  CODE l J )  =  I  CODE  I  IX  I 


MAIN  61 
MAIN  62 
MAIN  63 
MAIN  64 
MAIN  65 
MAIN  66 
MAIN  67 
MAIN  66 
MAIN  69 
MAIN  70 
MAIN  71 
MAIN  72 
MAIN  73 
MAIN  74 
MAIN  75 
MAIN  76 
MAIN  77 
MAIN  78 
MAIN  79 
MAIN  80 
MAIN  81 
MAIN  82 
MAIN  83 
MAIN  84 
MAIN  85 
MAIN  86 
MAIN  07 
MAIN  88 
MAIN  89 
MAIN  90 
MAIN  91 
MAIN  92 
MAIN  93 
MAIN  94 
MAIN  95 
MAIN  96 
MAIN  97 
MAIN  98 
MAIN  99 
MAIN  100 
MAIN1C1 
M  A  I  N  1  0  2 
M  A  |  N 1 0  3 
M  A  I  N 1 04 
M  A  |  N  1  0  5 
MAIN106 
M  A  I N 1 07 
H  A  |  N  1  0  8 
MAIN109 
M  A  I  N  1  1  0 
M  A  I  N  1  1  1 
M  A  I  N 1 1  2 

MAIM  13 
MA INI  16 
MAIN1 1 5 
M  A  I N 1 1  6 
MA  INI  17 
MA  IM  1  fl 
MAIN119 
MA  IN120 


F Yt J ) =F Y{ IX | 

MA  IN121 

MAIN122 
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CCNT i  \ue 

MA  1N123 

J  J* J  J ♦  1 

MAIN124 

KjR.  1 

H  A  I N 1  2  5 

CO  TO  20 

M.  A  I N  1  2  6 

22 

COM l NUE 

MA  INI  27 

IF  t FY (  1  I.CT . F Y ( 2  )  1  OC  TO  2  3 

MAIN128 

NFUNI  JJ  1  -F'M  2  1 

M  A  1  N  1  2  9 

1C0UE  (  J  J  !«.  ICOTE  (  ;  ) 

MA  INI  30 

Nf UN (  JJ  ♦  1  I  -V.N  (  1  I 

M  A  I N 1  3  1 

icooe  (  j  j  *n  =  i  c  one  ( ;  i 

MAIN132 

GO  TO  24 

M  A  I  N 1  3  3 

23 

f  FCNI  JJ  ]  -"Ml  1  1 

MAIM  34 

iccoe  (  J  j  i  Me  onr  1 1 1 

MAIN  35 

f;FuN  l  J  J  *  I  )  N  (  2  1 

M  A  I  N  1  3  6 

1  COOE  <  J  J  *  l  >  -  1  r  C.ur  (  2  ) 

M  A  I  N 1  3  7 

24 

CONT I  NUE 

MAIN139 

25 

NUPs  NFllN  (II) 

M  A  I  N 1  3  4 

fao=\f  uni  1 1  ♦  i  1 

M  A  ( M  4  0 

xi  {  I  )  =  xv.p(  1  ! 

MA  INI  4  1 

XU ( |  ) =XOP(  (  ♦  1  ) 

M  A  IN  1 42 

fUPiO. 

ft  I  N 1 4  3 

c 

M-  A  I  N 1  4  4 

c 

INTEGRATE  flOOY  SECTION  10  FIND  “ASS*  AM)  ACCUMULATE. 

M  A  I  N  1  4  5 

c 

M  A  I  M  4  A 

CAH  RMF(  IN  1  f  X,F  I  ,  X>.  I  I  i  >  XIJ(  |  1  .  1  0,  *  •  (  -61  ,  C.  ) 

MA|N147 

f  L  I  (•  =  1 . 0 

M  A  I N  1  4  B 

XM»XM<P  1  *R0(  fi'.iP  !  •  (  I 

M  A  !  N  1  4  9 

c 

M-A  INI  50 

c 

INTEGRATE  (UiOY  SlClION  TO  M NO  C.C,  10CATICK. 

M  t  IN  1  5  1 

c 

M  A  I  N  1  5  2 

CALL  RMIIC’IN  (  FX,  1  (  1  ,  XL  (  1  1,  XIJI  I  1  ,  1C.  *  *  (-6)  ,  C.  ) 

MAIM  53 

f L!P*2.0 

M  A  1  M  S  4 

XCGM  F  !  / f  I 

MAIM  5  5 

c 

MAIM  56 

c 

1  N  T  f  G  R  A  T  L  r.OOY  SfCTION  ’0  (IN!)  A  X  |  A  1.  MOM^M,  ANG  ACCUMULATE. 

MAIMS  7 

c 

M  A  I  M  5  3 

CALI  RMRCIN  ( F <* C 1 , XL (  1  )  .XU!  I  1  . 10. I-t  1  . C,  1 

M  A  I  ft  1  5  3 

FLIP, 3.0 

M  A  I Nl 6 0 

A  “I  .  A  “  I  *  0,  c.c  *>  1  o--  L  1  Ni;P  I  •  U  1 

M  A  I  N  I  6  1 

c 

M  A  l  e;  1  b  ? 

c 

INTEGRATE  ejOY  SfClIUN  Tu  M  NO  IKANSVlRSE  MOMENT,  AND  ACCUMULATE. 

“  A  |  N  1 6  3 

c 

M  a  I  N  1  4  4 

Call  Rff-GIN  (  1  x  ,  Of,  1  ,)  1  (  1  ) ,  V:  (  1  ) .  1  n.  <  *  ( -6)  ,  0.  ) 

M  A  I  N  L  (,  5 

2  f  =  /  1  *  f  |  •-  <t  .  » (’  l  *  <  u  i  -.t;p  1 

MA  INI  A'; 

K il  =  K  H  ♦  P  1  »l>0 1  Nl  (’(*(■)./  v  *0  I  1  ) 

M  MNU7 

IF  I  |  I.LU.  J  1  )  CO  TO  .»<> 

M  A  1  N  l  6  P. 

11=11*1 

M  A  I M  6  9 

GO  1 0  2  S 

MA IN  170 

26 

1  M  *  1 

V  A  I N  1  7  I 

If  (I.EC.L)  GO  TO  27 

MAIN172 

(.0  TO  IN 

ye  inn  7 

c 

MAIM  7  4 

c 

CAltUlAll  (..1,.  I'.CAIir’;  AM)  I.HAl  TKA.NSVTKSC  MOMENT. 

M  A  I  M  7  N 

t 

M  A  I  f;  1  7  6 

2  7 

(  Cl’-;  T;  J  -  /  t  ft". 

MAIM  7  7 

F.  !•))  -  XM»  f  MV.;- J  ♦  •  ; 

M  A  I  N  1  7  % 

c 

K  A  I  N  1  /  9 

PxIN’J  OiJl  M SUITS,  Af-,0  P.LU;"N  f  ■;.<  A  M  W  CASE. 

M  A  IMHO 

k 


% 


*£• 


I 

s 

\ 


MA  IN181 

WRITE  (6,30)  XM,CGPROJ,ANI,B,|ACODE(I),t»l,t)  NAIN192W 

GO  TO  l  MAIN183 

MINIS'. 

29  FORMAT  ( 6F10.6.2A9, 12)  MA1N185 

29  FORMAT  (/5X.5H  MASS,  1  2  X  »  2hCG*9X,6HAX  MCM,  IX  ,  9H  TRANS  MOM  ,  1  1  X  ,  AHCODEM.UN1 86 

1/)  MAIN107 

30  FORMAT  (A(2X,E12,5I  , 3X ,2 A9 I  MAIN188 

31  FORMAT  (////7X,lHA,12X,lMA,12X,lHX,l2X,lHV,l2X,lHR,12X,lHXf12X,lFlXNAIN189 

1 ,9X,7H0ENS1TY, 10X,enC0MM6NTS/8X, IHC,  l?x, KHl , 12 X , 1HC , 1 2X , l HC , 25 X, 1HMAIN190 
20,12',  |HF / )  M A  IN  19 1 

32  FORMAT  (1H  ,8613.5.3X,2A9)  MA1N192 

ENO  MAIN193- 

FUNCTION  FX  (X)  •  193*  1 

DIMENSION  XC(AO),  YC(AO),  RC(AO),  AOIAO),  AltAO).  XXCIACl,  Fx  2 

1  ICOCEtAO)  FX  3 

COMMON  XC,  YC.  RC,  AO,  Al,  NLP,  NLO,  FLIP,  ICOCE  FX  A 

XXC( NUP) =  X  — X  C ( NUP )  FX  5 

IF  (RC(NUP). EC,0.  I  XXC(NUP)»0.  FX  6 

XXC(NLO)=X-XC.(NLCi  FX  7 

IF  (RC(NLO). EG.O.  )  XXCINLOMO.  Fx  8 

IF  I  ICOOF  (NUPI»EC*-1I  GO  TO  1  FX  9 

YU«IYC(NUP)t (RC(NLP|**2-XXC(NUR)**21 »*0.5»A0(NUP)*Al(NUP) *X) «*2  FX  10 

GO  TO  2  FX  11 

1  YU«( YC(NUP>- I RC(NUP I »*2-XXC (NUP I *»2 ) **0.5* AO(NUP) *Al( NUP) *X ) **2  FX  12 

2  CONTINUE  FX  13 

IF  (  ICOOE (NLO).EC.-l )  GO  TO  3  FX  IA 

YL«(YC(NLOIHRC(NLOI**2-XXC(NLOI**2)**C.5*AO(NLOI»A1(NLO>*X)«*2  FX  15 

GO  TO  A  FX  16 

3  YLa(YC(NL0)-(RC(NL0)**2-XXCIM.0)A*2I**0,5*A0(M0I*A1(NLC)*X)**2  Fx  17 

A  CONTINUE  FX  18 

FX»YU-Yl  FX  19 

IF  1FUP.EQ.0.)  RETURN  FX  20 

IF  (FUP-2.)  5.6,7  FX  21 

5  FX«X*FX  FX  22 

RETURN  Fx  23 

6  FX*YU**2-YL* *2  FX  2A 

RETURN  FX  25 

7  FX=FX*(x*»2)  FX  26 

RETURN  FX  27 

END  FX  28- 

SUBROUTINE  RM8GIN  (  FX.Fl  ,Ll,U(.  ,TOL  ,PC  I  *  221*  2 

REAL  LL  RMBGN  2 

DIMENSION  A(  9  )  »  0(91  RM8GN  3 

00  1  1=1,9  RMBGN  A 

A (  I  1=0.  RMBGN  5 

1  6(11*0,  RMBGN  6 

XL  *L  L  RMBGN  7 

F A  =  F X ( XL )  RMBGN  9 

FaFX(UL)  RMBGN  9 

H«UL-XL  RMBGN10 

A(  1)=,5*H*(FA«F)  RMBGNll 

IP*1  RMBGN12 

IC=0  RM0GN13 

IS=1  RMBGN1A 

IF  I  PC, EO.O, I  GO  TO  2  RMBGN15 

WRITE  (6,11)  H,  !C,(A< | |,  1  =  1, Al  RM8GN16W 

2  IC=1  RMBGN 1 7 

3  Hl*H  RMBGN18 

H».5*H  RMBGN 19 


X  >XL  ftH 

RRBGN20 

suft=o. 

Rft8GN21 

DO  A  1  - 1 , IS 

RMBGN22 

SUM*FX( X ) ♦ SUM 

RRBGN23 

4 

X=Hl*X 

RB8GN24 

I  S  *  l  S  ♦  1  s 

RRSGN25 

B  t  l  ) *• 5* ( A  (  1  MHUSUH  I 

RM8GN26 

C«4. 

RRBGN27 

CO  5  J*  l  «  I  P 

RRBGN23 

K  *  J  ♦  l 

ft  M  8GN29 

S ( K )  =  ( C • 6  C J)-A(J)  1/IC-l.  I 

RftBGN30 

5 

04. *C 

RftBCN3l 

IF  ( PC.tC.O. I  GO  TO  6 

RftfiGN32 

WAITE  (6,111  H,  IC.  I  B(  II.  I  ■  1 .01 

RftOGN33w 

6 

CO  7  Jit , IP 

R  ft BGN34 

K  *  J  ♦  1 

RBBCN3S 

ABC  =  ABS (  (R(J)-H (Kll/UIKI  ) 

RRBCN36 

IF  (  ABC-TOL.  I  f  .!),  1  CO  fO  1  0 

ft  RBGN37 

ABO  AB5(  (  A  (  K  I-ni*.  »  »  /P.  (K  II 

ft  ft  B  G  N  3  8 

IF  (ABC-TCL.LE.C. I  CO  TO  10 

R  BPCN3  9 

7 

CONT  I  SUE 

ft  ft  PC,  N4  0 

IF  (  IP.cC.ei  GC  TO  8 

R  ft  0 GN4 1 

1  P  =  I  F»  ♦  1 

R  ft  BCN4  ? 

e 

I C  “  I C  ♦  1 

R  ft  BGN4  3 

CO  <9  JM.9 

R  ft  DGN4  4 

9 

A ( J  )  =B( J ) 

Rft  EGN45 

IF  <  1C.U.  10  1  GO  TO  3 

R  ft  BGN46 

WRITE  (6.121 

R  ft  0.  G  N4  7  w 

10 

f  i  =  e  (  k  i 

R  ft  BGN4  R 

RF.  TURN 

R  ft  2 GN4  1 

R  ft  6  0  N  6  0 

11 

FORMAT  1  1?E1 4. 7,  I  A , Of  12. 6 1 

RfthGNSl 

12 

FORMAT  I37H  RVRGIN  CIC  NOT  CONVERGE  IN  10  STEPS. 1 

ft  ft  RGN5? 

END 

R ► 6GN53 - 

SUBROUTINE  PKM  (X.N.XMAX, Jl 

•  274*  3 

DIMENSION  X  (  5  001 

P  ft  A  x  2 

X  ft  A  X  =  X  (  1  1 

P  ft  AX  3 

J  *  1 

Pft/X  4 

CO  2  l=2,N 

P  ft  A  X  5 

If  l XMAX-X I  lit  1,2,2 

P  ft  A  X  6 

l 

J=  1 

P  ft  AX  7 

XM AX  =  X(  1  ) 

p  ft  a  x  e 

2 

CONI INUE 

P  ft  A  X  9 

RE  TURN 

P  ft  A  X  10 

f  ND 

P  ft  AX  li¬ 

SUBROUTINE  P  M  |  S’  (X.N.XMIN.JI 

ft  285*  4 

DIMENSION  X ( 5001 

Pft  IN  2 

J  =  1 

Pft  IN  3 

X  M  1 N  *  X  1  1) 

Pft  IN  4 

00  2  1  =2  ,N 

PftIN  5 

IF  (  XMlN-Xl  11)  2,2.1 

Pft  IN  6 

i 

J=  I 

PftIN  7 

XM IN- X ( 1 1 

PftIN  8 

2 

CONT I NUE 

PftIN  9 

H  t  TURN 

PftIN  10 

END 

PftIN  11- 

ENO 

DATA 

t 


rJ*A 
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